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Abstract 

This paper presents a global thermal model of the energy conversion of the 10 kW e i Eurodish dish/Stirling unit erected at the CNRS- 
PROMES laboratory in Odeillo. Using optical measurements made by DLR, the losses by parabola reflectivity and spillage are calcu¬ 
lated. A nodal method is used to calculate the heat losses in the cavity by conduction, convection, reflection and thermal radiation. A 
thermodynamic analysis of a SOLO Stirling 161 engine is made. The Stirling engine is divided in 32 control-volumes and equations of 
ideal gas, mass and energy conservation are written for each control-volume. The differential equation system is resolved by an iterative 
method developed using Matlab™ programming environment. Temperature, mass, density of working gas, heat transfers and the 
mechanical power are calculated for one Stirling engine cycle of 40 ms and for a constant direct normal irradiation (DNI). The model 
gives consistent results correctly fitting with experimental measurements. 

© 2008 Elsevier Ltd. All rights reserved. 
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1. Introduction 

Since July 2004, a 10 kW e i dish/Stirling unit is in opera¬ 
tion at the CNRS-PROMES (French National Centre for 
Scientific Research) laboratory in Odeillo in France. This 
system is one of the several country reference units of the 
envirodish project. It is a Eurodish system developed by 
the deutsches zentrum fur luft und raumfahrt (DLR) and 
schlaich bergermann und partner (SBP) for solar electricity 
generation using a Stirling engine externally heated by con¬ 
centrated solar radiation (Keck et al., 2006). During these 
two years, the system has accumulated 2500 operation hours 
and 14.6 MWh el of electricity production with a power 
record of 11.1 kW e i at 974 W/m 2 direct normal insolation, 
which corresponds to a net solar-to-electric efficiency of 
21 . 6 %. 

The present work is focused on the thermal modeling. 
The objective is to evaluate the optical and thermal losses 
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of the system components, to identify the key parameters 
of the design and to study their influence on the perfor¬ 
mances of the components according to various operation 
conditions. The dynamic simulation of the system is also 
targeted in order to evaluate the efficiency on a daily, 
monthly or annual basis. 

The thermal model is composed of a radiation transfer 
model for the cavity, which is coupled to the solar flux dis¬ 
tribution, and on a thermodynamic model for the Stirling 
engine. The results are compared to experimental data 
issued from a campaign of optical and power measure¬ 
ments made by DLR in 2005 under a direct solar irradia¬ 
tion of 906 W/m 2 (Reinalter et al., 2006). 

2. Methodology 

The system is divided into five components, each of them 
representing a stage of the conversion cascade (Fig. 1). The 
optical losses of the concentrator (reflectivity, spillage) are 
determined using the distribution of solar flux density in 
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Nomenclature 

Aw heat transfer area, m" 

C p massic heat capacity at constant pressure, 

J kg -1 K -1 

C v volumic heat capacity at constant pressure, 

J kg -1 K -1 
D diameter, m 

Diss thermal dissipation, J 

DNI direct normal insolation, W m _ “ 
k thermal conductivity, W K -1 m _1 

/ receiver tilt angle, rad 

fr friction factor 

F view factor 

Gr Grashof number 

^ i 

h convective heat transfer coefficient, W m -_ K“ 

J radiosity, W m _ “ 

L length, m 

m working gas mass, kg 

m mass flow, kg s -1 

M molar mass, g mol -1 

n number of regenerator screens or absorber tubes 

Nu Nusselt number 

P pressure, Pa 

5Q heat power in Stirling engine, W 

R ideal gas constant, J kg -1 K -1 

Re Reynolds number 

R t thermal resistance, KW” 1 

S area, m 2 

S SQC section area, nr 

t time, s 

T temperature, K 

th thickness, m 

'j 

V volume, m 

v velocity, m s -1 

bW mechanical power, W 


Greek symbols 

n T Stirling cycle period, s 

s total hemispherical emissivity 

cp solar flux, W m _ “ 

G heat power in the receiver, W 

A P pressure drop, Pa 

(j) phase angle, rad 

<y a 

cj Stefan-Boltzman constant, W m _ “ K“ 

X wavelength, pm 

'j 

p m mass density, kg m“ 

p reflectivity 

co angular speed, rad s -1 

Subscripts 
a aperture 

am ambient air 

C compression control-volume 

cav cavity 

cer ceramic 

d dead 

inc inconel 

E “expansion” control-volume 

H “absorber” control-volume 

par parabola 

sec section 

w wall 

Exponents 
Cd conduction 

Cdv conduction and convection 

Cv convection 

LT low temperature, spectrum radiation emitted by 
a black body at a temperature about 1100 K 
rad thermal radiation 

S solar spectral band (X < 3 pm) 


planes close to the aperture (Reinalter et al., 2006). For a 
given direct solar irradiation, the thermal model of the cav¬ 
ity and of the absorber estimates the energy provided to the 
working gas of the Stirling engine and the heat losses by 
reflection, thermal radiation and convection. A nodal 
method is used to simulate one Stirling engine cycle, yield¬ 
ing the estimation of the generated mechanical power. 
The generator efficiency and the consumptions of all single 
components (water pump, tracking system, cooling fan) are 
user-defined inputs. 

3. Optical and power measurements 

The effective parabola area considering the shadows is 
52.9 m“ and the mirror reflectivity is measured at 92.5%, 


a value somewhat lower than the design value (94%). Rein¬ 
alter et al. (2006) assessed the optical performance of the 
concentrator under elevated direct solar irradiation 
(906 W/m 2 ), using the flux mapping system developed by 
DLR and described by Ulmer et al. (2002). The solar recei¬ 
ver is composed of a ceramic cylindrical cavity 30 cm in 
diameter and 12 cm in depth, the concentrated solar radia¬ 
tion enters through an aperture of 19 cm in diameter. The 
cavity walls are insulated using a silica-based ceramic mate¬ 
rial. The hexagonal absorber composed of 78 tubes made 
of inconel of 3 mm outer diameter is placed at the bottom 
of the cavity. 

Fig. 2 shows the solar flux map in the focal plane nor¬ 
malized to 1000 W/m 2 of direct normal irradiation (DNI) 
and 94% concentrator reflectivity. The aperture of the 
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solar energy 




cavity diameter 


solar flux 
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Fig. 3. Solar flux map in the absorber plane. 


cavity is located in the focal plane, it has a diameter of 
19 cm and it intercepts 85% of the incoming concentrated 
solar energy. 

Fig. 3 shows the solar flux distribution in the absorber 
plane, which is located 12 cm behind the aperture. It is 
observed that 78% of the concentrated solar energy directly 
hit the hexagonal absorber, while the remaining 7% hit the 
cavity wall. The flux distribution is rather inhomogeneous 
with a important peak of 1583 kW/m 2 and the mean solar 
flux on the absorber is measured at 702 kW/m 2 . 


4. Receiver model 


A nodal method is used to calculate the heat losses by 
reflection, by thermal radiation and by convection out of 
the cavity, and by conduction through the ceramic walls. 
The energy provided to the working gas flowing through 
the absorber tubes is also calculated over one complete 
cycle of the Stirling engine. 
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Fig. 2. Solar flux map in the focal plane. 


The receiver is divided into 11 control-volumes, named 
i = 1 to 11: 


- 8 for the absorber (i = 1 to 8) to take into account the 
inhomogeneous flux distribution, 

- 2 for the irradiated and shadowed ceramic cavity 
walls (/ = 9 and i = 10, respectively), 

- 1 for the aperture (i = 11). 

Each insolated control-volume (6 for the absorber, and 
1 for the ceramic walls) is assumed to receive a uniform 
solar flux density calculated from flux map in the absorber 
plane (Fig. 3) and given in the Table 1. The highest solar 
flux density of 1068 kW/m", a value lower than the flux 
peak of 1583 kW/m 2 is located on the “absorber” con¬ 
trol-volume 5. 

Fig. 4 shows a cross-section of the receiver and Fig. 5 the 
8 control volumes of the absorber. Note that both control- 
volumes 1 and 8 are not insolated; they are hidden behind 
front tubes. 

For each control-volume /, the energy balance is written 
with one equation (Eq. 1): 

(P ” C P =E 0 V-‘ +0 -' 

j 

where 0 y _i is the heat transferred by conduction, convec¬ 
tion or radiation from each control-volume j (j - 1 to 
11) to the control-volume i, and (9 ? is the concentrated solar 
energy hitting the control volume i. 

4.1. Radiation transfer model 

The radiation transfers between the control-volumes are 
calculated using a radiation balance applied to the solar 
receiver. The spectrum of electromagnetic radiation is 
divided in two bands in order to take into account the spec¬ 
tral dependence of the optical properties of inconel and 
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Table 1 

Measured mean solar flux for each “absorber” control-volume 


Control-volume 


1 

2 

3 

4 

5 

6 

7 

8 

Absorber 

Mean solar flux 

kW/m 2 

0 

501 

403 

1052 

1068 

569 

521 

0 

702 




Fig. 5. Geometry and control volumes of the Eurodish absorber. 


Table 2 


Inconel and ceramic properties 



8 

P 



1073K 

p s (X < 3 |am) 

p LT (X > 3 |am) 

Inconel 

0.88 

0.07 

0.14 

Ceramic 

0.9 

0.15 

0.12 


ture is considered as a black surface at ambient air 
temperature. The control-volume surfaces 1 and 8 are hid¬ 
den and they do not participate in the radiation balance. 
For each control-volume surface i, radiosities in the solar 
spectral band and in the LT spectral band are, respectively, 
given by Eqs. (2) and (3). Total radiosity J h given by Eq. 
(5) is the energy flux leaving the surface i composed of 
emitted plus reflected energy. In the solar spectral band 
the radiosity is only composed of a reflected solar energy 
and in the LT spectral band the radiosity is build from 
the emitted energy and the reflected LT energy by the other 
control-volume surface i. 



li 


J LT = s,(jTI + pl T 


W l 


E F ” J 


LT 

i n 




( 2 ) 

( 3 ) 


The view factors Fy are computed by integration of dif¬ 
ferential equations for discrete areas issued from Spiegel 
and Howell (1992). 

Knowing the radiosity, Eq. (4) gives the net radiation 
balance 0) ad for each control-volume: 

&f = S i {J i -cp i - jr FijJj) (4) 

7=1 lj '^8 


where J is the total radiosity of the surface /: 

J = J S +J LT (5) 

The net radiation balance of the aperture, 0™ , repre- 
sents the receiver heat losses by reflection and by thermal 
radiation. 


ceramic. The first one (S) corresponds to the solar radiation 
spectrum (2 < 3 pm) and the second one (LT) to the spec¬ 
trum radiation emitted by a black body at a temperature 
about 1100K (2 > 3 pm). The optical properties of inconel 
and ceramic are given in Table 2. 

For these two spectral bands, the control-volume sur¬ 
faces are assumed opaque diffuse and grey. The cavity aper- 


4.2. Conductive and convective heat losses in the cavity 

Some works in the literature provide correlations to cal¬ 
culate Nusselt number in a cavity as a function of the Gras- 
hof number Gr, receiver tilt angle /, aperture diameter Z> a 
and temperature walls T w . The natural convective heat 
transfer coefficient /? cav used in the receiver model is deri- 
vated from several correlations (McDonald, 1995). The 
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Stine and McDonald correlation for the Nusselt number 
Nu is given as follows: 


Nu = 0.088Gr 1/3 





where 


Rt“ = 


Lj + Lt 


1 


^inc^sec j—i 


and 



and 

s = 1.12-0.98 ( — ) (7) 

Vf-'cav ) 

The natural convective heat losses 0, v through the aper¬ 
ture are given by Eq. (8): 


S 


sec j—i 


pi_ 

4 



2 

outer 


D: -) • Yl tubes/ 


inner, 



4.4. Receiver heat losses 

Finally, the energy balance for the aperture control-vol- 
ume, 11, gives the total heat losses in the receiver. 


r -0 


flosses _ / 


rad 

11 


10 


+ E©;-n 

3 = 1 

10 

+ E © ; C :,1 

7=1 


cav 


Reflection and IR emission losses 

Conduction losses through the ceramic cavity walls 

Convection losses 



®f V = J2 h ™ S l T wi-T’am) 


( 8 ) 


i¥^ 11 


where /z cav is derivated from Nusselt number: 


(9)Acav = 


Nuk 


am 


L 


cav 


The conductive heat losses 0 Cdv through the ceramic 
receiver walls and then evacuated to the ambient air by nat¬ 
ural convection are calculated from thermal resistance 
Rt Cdw . jj ie ceram j c thermal conductivity k cer and the natu¬ 
ral convective heat transfer coefficient /z am are assumed 
constant. The expression of conduction losses is given by: 

(T wi - T am ) 


0 cdv = 

f^ cdv 

where for the control-volumes 9 and 10: 
j^cdv_ 1 _|_ 1 ^ / D cav + 2th cer 

h‘\ m S / kcrrlnL 


(9) 




cer' 


D 


( 10 ) 


cav 


and for the control-volumes 1 to 8: 


Rt Cdv = 


1 


+ 


th 


cer 


h'dm^i k CQY Si 


( 11 ) 


4.3. Heat transfers in the absorber 


The flux density of solar distribution on the absorber 
tubes is inhomogeneous inducing temperature gradients 
and low heat transfers by conduction. Conductive heat 
transfers 0^ along the thin of the tubes between the eight 
“absorber” control-volumes are calculated by Eq. (11). 


T — T 

/-vCd _ 1 wy 1 wi 


(12) 


4.5. Absorber energy balance 

For the absorber control-volumes, Eq. (1) can be written 
as follows: 


(Pm^p V) Hl 


d T 


wH/ 


j^i 


d t 


= - E 0 “ - 0 “>i - 0 f-n - S Qm 

j 

(16) 


where 3Q n is the energy provided to the Stirling engine, 
depending on the pressure, the temperature, the mass den¬ 
sity and the velocity of the working gas. 

5. Stirling engine thermal model 

A Stirling cycle machine operates on a closed regenera¬ 
tive thermodynamic cycle using a working gas, and subjects 
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Fig. 6. Control-volumes in Stirling engine. 
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the gas to expansion and compression processes at different 
temperatures. Since its invention by Robert Stirling in 1816 
different engine configurations have been developed 
(Thombare and Verma, 2008). The Eurodish Stirling 
engine, SOLO 161 is of the alpha-type. Alpha-type engine 
have two pistons in separated cylinders arranged in “V”, 
which are connected in series by three exchangers: a heater 
(here a solar absorber), a regenerator and a cooler. 

Many works presents or uses thermal model more or less 
complex to approximate the real behaviour of the alpha- 
type Stirling engine. A review of these thermal models 
can be found in Dyson et al. (2004) and in Ercan Ataer 
and Karabulut (2005). 

For this study, a nodal analysis of the SOLO 161 is per¬ 
formed. The Stirling engine is divided in 32 control vol¬ 
umes as shown in Fig. 6. The control volumes of 
compression V c and expansion V E spaces are variable 
and expressed by Eqs. (16) and (17). The other control-vol¬ 
umes have fixed volumes. 

Fe(*)=- y(sin(©?) + l) + F E ,d (17) 

Vc(t) =^Y(sin(cot- (f>) + 1) + V c ,d (18) 


In the control-volumes of the compression and expan¬ 
sion spaces, the mechanical power bW\ is expressed by: 

dV 

bWt = — P — y— (22) 

d t 

And in the “heat exchanger” control volumes, the heat 
transfer 5Q t is given by: 

SQt = hiAwi(T wi - T t ) (23) 

h[ is the forced convective heat transfer coefficient calcu¬ 
lated with two correlations. The classical Colburn correla¬ 
tion is used in the 8 “absorber” control-volumes and in the 
10 “cooler” control-volumes. For the 10 “regenerator” 
control-volumes, the correlation given by Lemrani (1995) 
is used: 


Nut = 0.42Re°' 56 

L l 



Dissi i is proportional to the pressure drop A P t depend¬ 
ing to the friction factor. As for the forced convective heat 
transfer coefficient, two correlations are used: 


PM 

Diss i = A Pi - 

RTi 


th \n hi out 



where: 




A P,{ 


\ 


friiP„,3 n i 


= ^ + 0.337 For the “regenerator ” control volumes 

r 1 P+mvrLi 

J ri 2 A 

I 

r_ 64 

Ret 

If Re, < 2000 

fn{ 

1 

For the “absorber ” and “cooler ” control volumes 



[ = 0.316Re“ 025 

If Re, > 2000 



Mass equation (18), energy equation (19) and state equa¬ 
tion (20) are written for each control-volume. 

The following assumption are made in order to simplify 
the model equations: 

(1) The working gas is assumed to be an ideal gas, 

(2) The total mass of the gas is constant within the Stir¬ 
ling engine, 

(3) The instantaneous pressure is uniform, 

(4) The compression and expansion spaces are adiabatic, 

(5) The gas flow is one dimensional and quasi steady. 


dnij 

dt 



til out) i 


bQi T Dissi H - T m m m Cp{j ou t ) T outtiiout)i T bJV \ 





6. Numerical resolution 

The differential equation system is resolved by an itera¬ 
tive method developed in Matlab. Temperature, mass, den¬ 
sity of working gas, heat transfers and the mechanical 
power are calculated for one Stirling engine cycle of 
40 ms and for a constant direct normal irradiation. Peri¬ 
odic steady cyclic condition is obtained after few iteration 
cycles. 

Mean value per cycle of the energy powers G are the 
integration of the discrete powers bG over the time and 
are calculated by Eq. (26): 

f . bGdt 

cycle 

“ nf 

Few input data are necessary: 

(1) Geometric data of the receiver and Stirling engine. 

(2) Working gas characteristics of the Stirling engine. 
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(3) Eurodish operation data: engine speed, DNI, total 
working gas mass. 

It is important to point out that Eurodish system regu¬ 
lates the total working gas mass and consequently the cycle 
mean pressure, to keep constant the maximal absorber 
temperature at 1053 K. In the Stirling engine model, the 
total working gas mass is an input data and the absorber 
temperature is an output data. To solve this problem, a 
loop on the total working gas mass has been added in 
the numerical resolution in order to obtain a maximal 
absorber temperature of 1053 K over one Stirling cycle. 

7. Results and discussion 

Our numerical study is compared to experimental mea¬ 
surement given by Reinalter et al. (2006). The Eurodish 
operation data during the experimental measurements are 
provided in Table 3. Fig. 7 gives the model results and 
the experimental measurements of the solar energy dis¬ 
patching. The “numerical errors” designation is the differ¬ 
ence between the energy input and output of the Eurodish 
system: For the experimental measurements this term is 
due to measurement errors (Reinalter et al., 2006) and 
for the model results, it includes errors coming from 
numerical integrations. Fig. 7 shows the considerable 
energy losses in the receiver in particular by spillage. The 
losses by spillage, conduction, convection and radiation 
are estimated at 13.5 kW, which represents 28% of the solar 

Table 3 


Eurodish operation data 


Working gas 

Hydrogen 


DNI 

W/m 2 

906 

c 

‘-’par 

m 2 

52.9 

Solar energy 

kW 

48 

P 

bar 

[130-140] 

T wH max 

K 

1053 

T e 

K 

[903-923] 

[T n —T out ] water cooler 

K 

295-305 

Tam 

K 

293 

Nt 

rt min -1 

1500 



Fig. 8. Absorber temperature distribution: model results in underlined. 



Fig. 9. Power provided to Stirling engine, power provided to the cooler 
and electrical power as function of DNI. 

energy incoming of the dish; 31 kW is provided to hydro¬ 
gen and 19 kW is evacuated to the cooler. The developed 
model gives good results in agreement with the measure¬ 
ments particularly in the receiver. The efficiency of the 
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□ Cooler 
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Fig. 7. Solar energy dispatching, left: model results, right: experimental measurements (Reinalter et al., 2006). 
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0.5 7.5 



39.0 



□ Dsh reflection 
II Spillage 

□ Conduction 

□ Convection 

■ IR+Reflecticn 3.0 

□ Cooler 

□ Friction dissipation 

□ Electrical dissipation 

■ Consumption 

■ Electricity 

□ Numerical error 


13.9 


33.9 


10.9 


11.0 


Fig. 10. Solar energy dispatching for two DNI values, left: 800 W/m 2 , right: 400 W/m 2 . 


receiver defined as the proportion of the thermal energy 
introduced into the engine to the power that enters in the 
cavity is estimated at 81.4% and is measured at 82.6%. 
The observed differences between the model results and 
the measurements are due to the Stirling model. The Stir¬ 
ling engine efficiency, which represents the conversion of 
thermal to mechanical power is calculated at 34.4% and 
is measured at 39.2%. This difference can be explained by 
correlations used to calculate the forced convective heat 
transfer coefficients in particular in the regenerator where 
the gas flow is complex (Bonnet, 2005). Finally, the model 
finds a net solar-to-electric efficiency of 19.1% instead of 
22.5%, which gives a difference of 1.3 kW e i 

Fig. 8 shows the absorber temperature distribution cal¬ 
culated by the model and measured by 20 thermocouples 
placed behind the absorber. The total hydrogen mass 
within the Stirling engine is adapted to obtain a maximal 
temperature about 1053 K. This maximal temperature is 
located on the surfaces 4 and 5. Fig. 8 shows a good agree¬ 
ment between the model results and the measurements. The 
difference between the maximum and minimum tempera¬ 
ture is calculated at 130 K when measurements gives about 
100K. The temperature of the expansion space is estimated 
at 914K and the cycle pressure at 115 Bar. 

Fig. 9 shows the power provided to the Stirling engine, 
the power provided to the cooler and the electrical power 
as a function of the DNI. The maximal absorber tempera¬ 
ture is constant at 1053 K. The model results are compared 
to experimental values of DNI, net electrical power and 
temperature in the cooler made automatically every second 
on the Eurodish system (Reinalter et al., 2006). The exper¬ 
imental net solar to electricity efficiency increases as a lin- 
ear function of a DNI up to 800 W/m". When the 
insolation exceeds this level a blower in the receiver cavity 
is activated to reject additional heat and maintain the 
absorber at a fixed temperature of 1053 K. Model give 
good results. For the electrical power the difference with 
experimental values is less than 1 kW e i and is maximal at 
700 W/m". At 800 W/m" the power provided to the cooler 
is evaluated at 16.5 kW and is measured at 17.5 kW. This 
energy could be recuperated for a cogeneration applica¬ 
tion. In this case, the overall efficiency (Onovwiona and 
Ugursal, 2006) increases to 58%. For the power provided 


to the Stirling engine It can be seen a constant difference 
of 1 kW between the model and the experimental values. 

Fig. 10 gives the energy dispatching calculated by the 
model for two DNI values: 400 W/m 2 and 800 W/m 2 . 
The concentrator efficiency, defined as the proportion of 
energy entering in the cavity to solar energy intercepted 
by the dish is constant at 78.6%. The losses by conduction, 
convection and reflection and thermal radiation in the 
receiver increase from 5.2 kW at 400 W/m 2 to 6.5 kW at 
800 W/m 2 DNI, which gives a receiver efficiency of 72% 
and 80%. The Stirling efficiency is constant at 34.4% and 
finally the net solar electricity efficiency drop from 800 to 
400 W/m" due to the receiver efficiency decrease is 16.2%. 

8. Conclusion and perspectives 

This paper presented a thermal model of energy conver¬ 
sion of Eurodish dish/Stirling. The model compared to 
experimental measurements gives good results and permits 
to study the detailed heat losses. Currently, the net conver¬ 
sion solar-electricity efficiency is about 21% at a DNI value 
of 900 W/m 2 and at an ambient temperature of 20 °C. The 
model shows important heat losses in the cavity by spillage 
and radiation (reflection and IR-emission) and conse¬ 
quently an improvement potential of system efficiency. 

Future works will focus on the study of the influence of 
cavity physical parameters in order to design a new solar 
receiver. Using the ray-tracing code Soltrace and parabola 
slope error data, a parabola and receiver simulation are 
made to estimate and study the solar flux distribution in 
different planes. Then, theses flux maps will be introduced 
in the presented model. Furthermore, the model is used to 
evaluate energy provided to the cooler in order to test a 
Eurodish dish/Stirling for a cogeneration application. 
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